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It is demonstrated that thermal radiation of small occupation number is strongly nonclassical. 
This includes most forms of naturally occurring radiation. Nonclassicality can be observed as a 
negative weak value of a positive observable. It is related to negative values of the Margenau-Hill 
quasi-probability distribution. 
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The very signature of quantum mechanics, Planck's 
constant ?i, was measured for the first time in 1900 in 
blackbody radiation [1]. In order to derive the spec- 
trum for blackbody radiation, Planck had to assume that 
the energy of the radiation field could only assume dis- 
crete values. This was the first indication that the clas- 
sical electromagnetic theory of radiation did not provide 
a fully satisfactory model of nature. However, Planck 
himself maintained that the radiation field could be de- 
scribed by classical electromagnetism. For him, it was 
the emission from the material oscillators that was quan- 
tized [2] . More recent work has confirmed that the black- 
body spectrum can be explained in terms of stochastic 
electromagnetism [3]. 

The invention of the laser around 1960 sparked the de- 
velopment of a whole new are of physics, quantum optics. 
It was found that coherent states, a non-thermal state of 
light, was a natural building block in a quantum theory 
of radiation [4]. Glauber and Sudarshan demonstrated 
that any density matrix could be expanded diagonally in 
terms of coherent states [4, 5], 
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The weight function P{a) is known as the Glauber- 
Sudarshan P-distribution. Glauber also proposed that 
nonclassical states, i.e., states that cannot be modelled by 
a classical theory, are those for which the P-distribution 
fails to be a probability distribution [4, 6]. Put differ- 
ently, any state which cannot be expressed as a classical 
mixture of coherent states is nonclassical. 

Among pure states, the coherent states are the only 
ones satisfying the Glauber classicality criterion [7]. 
However, it was recently shown that coherent states may 
display nonclassical properties in weak measurements 
[8]. This demonstrates that Glauber's classicality cri- 
terion has not taken into account the possibility of weak 
measurements. Weak measurements were discovered by 
Aharonov, Albert and Vaidman (AAV) in 1988 [9]. 

Thermal radiation is almost the only form of naturally 
occurring radiation, emanating from any object of finite 
temperature. It can be seen in such diverse areas as the 
cosmic background radiation, the light from stars and 
from earthly sources such as a fire or a light bulb. The 



P-distribution for a thermal state is 
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where (n) is the expected photon number. This distribu- 
tion is neither negative nor singular, hence this is essen- 
tially a classical state according to the Glauber criterion. 

In this Letter, we will demonstrate that also thermal 
radiation may display nonclassical properties in weak 
measurements. The term "nonclassical" means, in this 
respect, that the effect cannot be reproduced by a model 
where the observables are stochastic c-numbers. We will 
also show that the effect is related to negativity of the 
Margenau-Hill quasi-probability distribution [10]. 

In a standard, projective measurement, the uncer- 
tainty of the pointer is initially very small so that the 
pointer can distinguish different eigenvalues of the ob- 
servable [11]. In the weak measurement scheme of AAV, 
the pointer is assumed to be in a pure, gaussian state of 
large uncertainty. In this way, the pointer cannot discern 
the different eigenvalues of the observable. The theory 
of weak measurements has been generalized to arbitrary 
states of object and pointer in Ref. [12]. Generally, a 
weak measurement of an observable c conditioned on the 
outcome of projective measurement of an observable d 
yields the real part of the weak value 
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Here p is the density operator of the object, c^, is in- 
dependent of the specific choice of pointer state. The 
pointer may be in an arbitrary state, pure or mixed, 
provided that the current density vanishes. Also, it is 
required that the interaction between the object and 
pointer should be sufficiently weak [12]. 

The Hamiltonian of the free radiation field is defined 



as 
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There are two contributions to the Hamiltonian from 
each of the quadratures q and p. In this sense, the energy 
of the radiation field consists of one energy contribution 
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from each quadrature. For a material oscillator, one term 
is kinetic and the other is potential energy. We shall con- 
sider weak measurements of the energy contribution from 
one quadrature postselected on the other. It was recently 
shown [8] that the weak value of postselected on g is a 
conditional moment of the standard ordered distribution 
[13] 
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where the standard ordered distribution is defined as [13] 
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The standard ordered distribution is the complex conju- 
gate of the Kirkwood distribution [14]. 

By using Eq. (1), we may express the standard ordered 
distribution in terms of the P-distribution, 
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By using the quadrature representation of a coherent 
state I a) [15] 

-■^ + V2aq- - 



\ -1/4 

a) = TT ' cxp 



1 



' 2 



and by inserting the P-distribution for a thermal state 
in Eq. (7), it can be shown that the standard ordered 
distribution for a thermal state is 
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where 
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is the variance of each quadrature. 

Thermal radiation is often well described by the spec- 
tral distribution of blackbody radiation. For blackbody 
radiation the expected occupation number is 

(") = — rij^- (11) 

At maximal irradiance of the Planck distribution, as 
given by Wien's displacement law, the occupation num- 
ber is of the order 10~^. The Margenau-Hill distribution 
[10], which is the real part of the standard ordered dis- 
tribution, has been plotted in Fig. 1 for an occupation 
number of 10^^, and in Fig. 2 at an occupation number 
of 1. 

The standard ordered distribution yields correct 
marginals when integrated over. Thus, the marginal dis- 
tribution for the quadrature q is 

{q\p\q)^ I dp S{q,p) = -L- e-'V(2.^). (12) 




FIG. 1: The Margenau-Hill distribution for thermal radiation 
at maximum irradiance given by Wien's displacement law, 
i.e., for an occupation number of (n) « 10"'^. The distribu- 
tion is virtually identical to that of vacuum [8, 16]. There are 
distinctly negative regions, which lead to nonclassical behav- 
ior. 
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FIG. 2: The Margenau-Hill distribution for thermal radiation 
for an occupation number of {h} = 1. The negative regions 
have practically disappeared. 



This is a standard gaussian distribution with variance cr^ . 
By inserting Eq. (9) in Eq. (5) we find that the weak 



value of p is 
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This is an inverted parabolic curve which is negative for 
\q\ > Vc^ + 4cr6. Thus, for sufficiently large q, the weak 
value of p^ is negative for any thermal state. A classical 
radiation model cannot reproduce this result. The theory 
of weak measurements can also be extended to classical 
theories. It can be shown that if observables are treated 
as classical c-numbers, the weak value of an observable 
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FIG. 3: The probability of observing a negative weak value 

of for thermal radiation as a function of the occupation 
number (n) . At maximum irradiance of blackbody radiation, 
given by Wicn's displacement law, the occupation number 
is of the order 10~^, which gives a significant probability of 
nonclassicality. 



is just the conditional expectation value over a positive 
joint probability distribution [12]. Therefore, in a clas- 
sical model, the weak value of a positive observable is 
always positive. 

The probability of observing a negative weak value is 
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The result is 
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This function has been plotted as a function of (h) in 
Fig. 3. We see that the result is nonclassical essentially 
for (n) < 1. For occupation numbers typical of ther- 
mal radiation states there is a significant probability for 
observing a nonclassical negative weak value of p^. 

In principle, weak measurements employ the same in- 
teraction Hamiltonians as strong (or projective) measure- 
ments (to this end, compare Refs. [12] and [17]). If a 
strong measurement can be performed with a specific in- 
teraction Hamiltonian, then a weak measurement can be 
performed with the same interaction type provided that 
the interaction strength is sufficiently weak. 

It could also be possible to observe the nonclassical ef- 
fect found in this Letter for a massive particle in thermal 
equilibrium in a harmonic oscillator potential. If a detec- 
tor is placed in a distance exceeding \/a^"+4a^ from it's 
equilibrium position, a weak measurement of kinetic en- 
ergy should yield a negative value on average. However, 
this requires observation in the low temperature regime. 
For a temperature of 1 /iK, the required oscillator fre- 
quency is of the order 100 kHz. 

From Eq. (4) for the Hamiltonian, it can be shown 
that 



This suggests an alternative strategy for measuring the 
weak value of . It could also be done by performing a 
weak measurement of energy postselected on the quadra- 
ture q, and subsequently subtracting the squared quadra- 
ture [18]. 

A comment can be made on the usage of pointer states. 
It has recently been shown that weak measurements can 
be performed with arbitrary pointer states, with the only 
restriction that the current density of the pointer state 
must vanish [12]. This means that a weak interaction 
between a thermal pointer state and another thermal ob- 
ject state may yield a pointer reading which cannot be 
explained in terms of a classical c-number model of radi- 
ation! 

In conclusion, we have demonstrated that thermal ra- 
diation displays strong nonclassicality in weak measure- 
ments. The effect is present for low occupation num- 
bers, as found in most naturally occurring forms of radi- 
ation. A weak measurement of the squared quadrature 
postselected on the second quadrature may yield a neg- 
ative weak value. This effect cannot be reproduced by 
a classical model of radiation where observables are c- 
numbers. The nonclassical effect is related to negativity 
of the Margenau-Hill distribution for the thermal state. 
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